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Abstract 

We derive new integral representations for objects arising in the classical theory of ellip- 
tic functions: the Eisenstein series Ea, and Weierstrass' p and ( functions. The derivations 
proceed from the Laplace-Mellin transformation for multipoles, and an elementary lemma 
on the summation of 2D geometric series. In addition, we present new results concerning 
the analytic continuation of the Eisenstein series as _an entire function in s, and the value 
of the conditionally convergent series, denoted by E2 below, as a function of summation 
over increasingly large rectangles with arbitrary fixed aspect ratio'^. 

1 Introduction 

In this paper we revisit the classical theory of elliptic functions as developed by Eisenstein 
and Weierstrass. Both of these researchers represented the meromorphic functions appearing 
in their theories as summations over a given lattice of elementary pole functions of a pre- 
scribed order. Our fundamental observation is that pole functions may be represented by 
exponentially-damped, oscillatory integrals. These representations depend on the complex 
half-planes in which the singularities lie, and are natural variants of the classical Mellin, or 
Laplace-Mellin, formulas which are valid for isolated poles lying in the right half plane, where 
?R.{lo) > 0. In more recent times such integral representations have resurfaced in the develop- 
ment of fast multipole methods for the solution of Poisson problems; in this context they often 
are referred to as "plane- wave" representations ([6] and, more recently, [3]). A key feature of 
these integral representations is that the pole centers appear in the exponents of the integrands. 
As a consequence the lattice summations are transformed into geometric series which may be 
summed explicitly underneath the integral. The result is a new class of integral representations 
for the Eisenstein series and other meromorphic functions of Weierstrass' theory. 

A brief summary of the paper follows. In the first section we review the definitions of the 
Eisenstein series En and the Weierstrass^ functions p and ^. We will analyze a generalization of 
Eisenstein's series which we denote by E's, the differences being: first, we consider s = a + it (z 
C, and second, we define Eg as a limit over lattice squares of increasing size, a significant point 
when 3?(s) < 2 and the sums are not absolutely convergent. In addition, in this preliminary 
section we provide elementary derivations of the requisite plane-wave formulas for general pole 
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functions of the form /(w) = and a summation identity for a two-dimensional geometric 
series. 

In the next section we employ the plane-wave representations and the summation identity 
to derive an integral representation for Eg for the case 5R(,s) > 2. Integral representations 
for Eisenstein's i5„^naturally follow for s = n > 3. Subsequently, we derive an alternative 
representation for Eg as a contour integral from which we deduce that the sums Eg, defined 
unambiguously for 5ft(s) >2, admit an analytic continuation as an entire function to the whole 
of the complex plane. As a corollary, we prove the existence of a finite limit for £2- We 
discuss E2 and its relation to Eisenstein's, E2- As the summation processes defining these two 
conditionally convergent series are distinct, one expects different limiting values. We derive a 
closed form correction term for this difference. (Note that both Eisenstein's convention and 
ours give £^1 = £^1 = 0.) In the following section we derive analogous integral formulas for 
Weierstrass' p and ( functions. We conclude the paper with a brief discussion of these integral 
representations in relation to previous research in the theories of lattice sums, and elliptic 
functions. _ 

In this last regard we take a moment to mention here that our formulas for Es are the 
natural lattice analogues to the well-known representation for Riemann's zeta function (which 
we denote with the subscript so as to distinguish it from Weierstrass' function of the same 
name) 



For example, in the case of a square lattice we derive the following integral expression for the 
classical Eisenstein series 



where /c is a positive integer divisible by four {Eh{i) — otherwise). The similarity between 
(1) and (2) is clear. For more general lattices, we replace i by r, fc € N by s € C, and the 
single trigonometric ratio in (2) by a sum of analogous ratios denoted by /i(t, A) and /2(t, A) 
defined in (17) and (19). The general expression is given in theorem 4. 

We note that a subset of the results presented below appeared previously in a slightly 
different form [7]. 

2 Preliminaries 

In this section we review the definitions of the Eisenstein series and the Weierstrass p and C 
functions. Furthermore, we derive elementary lemmas concerning plane-wave representations 
and a geometric series identity, both of which we will use repeatedly in the subsequent sections. 
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2.1 The Eisenstein series and Elliptic functions 

We are given a general lattice A c C which we describe by generators /x, u: 

A = {m ■ ii + n- y\m, n e Z} 

where ^,,y are complex numbers such that the lattice ratio, r = y/n, is not real. The fun- 
damental parallelogram is the set, Aq = {a + /3t||q:| < 1/2, < 1/2}. Eisenstein began his 
investigations of doubly periodic, meromorphic functions with periods jj and u through the 
study his eponymous series (see, for example, [12] and [9]) 

N M ^ 

Ek = lim V hm V k > 1; (3) 

n=-N m=-M ^ ' ' 

the elimination of the term m = n = is implicit here and below. For the magnitude of each 
term in the summand we have the bounds 

Kp-, v) ^ 1 ^ Kjfi, y) 



s/m? + n2 \m-ii + n-v\ ^/rv? + n 

Applying (4) and elementary estimates one verifies that the series (3) are absolutely convergent 
for n > 3, and thus are well-defined functions of the lattice En = En{A). The same estimates 
indicate that the series En are absolutely divergent for n = 1 or 2, hence the limiting operation 
specified in (3) plays a non-trivial role in the definition of these sums. Eisenstein proved that 
the limiting procedure (3) yields finite values of E^ even in these cases. We observe that, in 
particular, Ei = 0. This is in keeping with the fact that, using a symmetry argument in the 
absolutely convergent case, one may prove that i?„ = for all odd n > 3. Therefore, from the 
point of view of convergence, the only "interesting" sum is £'2- 

Eisenstein was cognizant that the value of E2 depends on the choice of limiting procedures; 
and he derived many identities which connect his summation process for E2 to others ([12]). 
We choose yet a different summation convention and define Eg as the limit of partial sums over 
"lattice-squares" of increasing size. We generalize further in considering complex exponents. 
Specifically, we define Eg by 

Es = lim y ^, (5) 

\m\,\n\<K \ ' / 

which we consider, initially, for s = l,s = 2 and §?(s) > 2. For non-integer s we consider 
the branch of the function with a cut along the "negative diagonal" of the lattice, {z = 
—t{ii + v), t > 0}. For all z in the closure of this cut plane we have 

< arg(z) < 9 + 2iT, where , . 

e = Avg{-^i-iy). ' 

We further enforce the convention that points of the lattice lying along the diagonal are con- 
sidered symmetrically, 

1 1/11 



3 



Wc will return to this point later. Finally, wc note that the large K limiting convention defined 
in (5) is relevant only in the cases s = 1 or s = 2. 

We will derive integral representations for Eg. Naturally, restricting s to the positive 
integers, our formula yields an integral representation for the classical Eiscnstcin series En — 
En, for n > 3. As for the conditionally convergent series, it is straightforward to verify that 
El = directly from (5). However, for s = 2, it is not a priori obvious that (5) will be finite 
for arbitrary A(/i, i/). Furthermore, if finite, the relationship between its value and Eisenstein's 
E2 is not clear. The finite existence of the limit (5) will follow as a consequence of the integral 
representations for Eg,^{s) > 2. In addition, we derive a formula which connects our limiting 
value to Eisenstein's. Even more, we prove that Eg admits an analytic continuation to s e C 
as an entire function. 

Some fifteen years on the heels of Eisenstein, in 1862 Weierstrass commenced his study of 
doubly-periodic functions. For Weierstrass, the fundamental object was his p function which 
he defined as 

1 NT- / 1 1 



p{x,A) 



+ 



E 



(7) 



weA\{o} 

Standard estimates demonstrate that p is an analytic function of x, and well-defined function 

of A. From (7) wc see that the p function has a double pole at the origin and every lattice 
translate. Weierstrass defined his ( function as an indefinite integral of p and developed the 
following summation representation 



1 



C{x, A) = - / p(s, A)ds = ^ + 



{X- 
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(8) 



Again, the correction terms inside the summation balance the asymptotics of the translated 
poles so that the sum converges absolutely for a; in a compact set containing no lattice points. 
As C is meromorphic and has only a single pole within a fundamental parallelogram, it can not 
be doubly periodic (see, for example, [1]). Nevertheless, Weierstrass derived many properties 
satisfied by this function and its relatives. The absolutely convergent sums (7) and (8) will 
serve as the starting points for the derivation of the integral representations for p and ( below. 

We conclude this section with a brief comment on freedom in the choice of generators for 
A. Note that Es,p,( satisfy simple rescalings with respect to /x 



Eg{n,iy) 
pix\iJ,,u) 
C{x\fJ,,u) 



(9) 



where t = vf/j.. 

Furthermore, one observes that two pairs of generators {fx, v) and (/x', u') will give rise to 
the same lattice if and only if they satisfy a linear system of the form 



a h 
c d 
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with a,b,c,d G X,ad — be = ±1. The set of aU two- by- two matriees with integral entries and 
determinant plus or minus one is known as the unimodular group in two variables. In this case 
the lattice ratios are transformed by the fractional linear transformation 

, a + br 
c + dr^ 

known as a unimodular substitution. One may prove that up to rescaling and unimodular 
substitution, any lattice ratio may be represented by a unique r chosen from the following 
fundamental region [1] 

' -I < 5i(r) < i 

< > 0, (10) 

|t| > 1, 

^ if |t| = 1, then ^{t) > 0. 

In summary, without loss of generality, we restrict our analysis to the "inhomogeneous" 
functions, which are obtained from (5), (7) and (8) by fixing fj, = l,y = t, and consider 
A = A(t) with T satisfying (10). For convenience we omit the variables fi, v below and write, 
for example, p = p{z,t). 



2.2 Plane- wave representations and a 2D Geometric Series 

To facihtatc our derivations we define the truncated lattice A^ = {cjm_„ = m + nr] \m\, \n\ < 
K} \ {0}. We further group lattice points into four overlapping "quadrants" 

A^ = A^_.)UA^,+)UAf_,.)UA^__) 

defined by 

^U-,,)=Wm,n\ l<m< K, |n|<TO} 

A^-H) ={'^™,"l l<n<K,\m\<n} 

AK_.) = {w„^,n|-^<m<-l, \n\<-m} ^ ' 

_) = {uJm,n\ -K < n < -1, |m| < -n}. 

We recall from the discussion following (5) that for non-integer s the shared boundary between 
A(^ ,) and A^ ^ is identical to the branch cut (see figure 5). 
We have the following elementary lemma 

Lemma 1. Assume a complex lattice A(r) and the quadrants defined as in (11). An isolated 
singularity of complex order s, 3?(s) >0 with branch cut defined as in (6) may be represented by 
the following plane-wave integrals, each of which is valid in the appropriate quadrant determined 
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by the location of the point uj. 



1 1'°° 

^ — {■KS poo 

'nTs/2 poo 
s/2 poo 

s) Jo ^ ' ' 



Proof: 

We have the representation of the T function 

r(s) = / A"-^e-^dA, U{s) > 0. 
Jo 

Assume u) € iR+, in particular, uj — it, t > 0. As r satisfies (10) we observe that — tt < 6* < 
—tt/2 hence arg(a)'*) = 7r,s/2. Keeping this in mind, we rescale the integration variable by t, 
factor the —1 in the exponential, and multiply and divide by exp(i7rs/2) to obtain 

POO 

r(s) = iM A"-^e-^*dA 
Jo 

= e-''^^/^^" / A^-ie'^'^dA. 
Jo 

Dividing both sides by r(s)a;* gives the desired result for uj = it. In a similar manner we 
prove the formula for u lying on any of the principal coordinate rays emanating from the 
origin, lo G ±R+,±iR+. The full expressions (12) then follow by analytic continuation into 
the appropriate quadrants, f 



We note that for integer s, the integral expressions may be continued further and are valid in 
the appropriate half-planes ±5f(a;) > and ±9(a;) > 0. 

Next we turn to our summation convention (5). From lemma 1, it is apparent that no 
single plane- wave expansion formula will be valid for all terms in the summands (5), (7), and 
(8); terms miist be grouped with respect to quadrant. As with the our convention of splitting 
contributions from points LO-m,-m lying on the cut in (5) equally between branches, we wish 
to treat each quadrant as symmetrically as possible. We define the symbol 6^71 for m,n G Z 
by ^ 

m = ±n 



2' 

1, otherwise 



(13) 
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By convention, we sum over the terms in the ^^-quadrant as 

K m 

J2 -^('^) ^ J2 J2 emr^/(w„^,„) (14) 

/.i(=A^ m—ln— — m 

K m-1 ^ \ 

= X] I 2^i^rn-m) + fi'^rn,n) + ^f{l^m,m) 1 (15) 

m=l \ n——m+l / 

We recognize the slight abuse of notation in (14) in that the left hand side is not interpreted as 
a standard sum over the set , but rather a modified sum in which the "diagonal" terms 

are added with a factor of 1/2; this is made explicit in (15). We choose this abuse of notation 
so as to not burden our summation symbols with a cascade of modifiers, and anticipate that 
it will not cause confusion below. 

The sum over the lattice square is the sum of the quadrant sums as in (14); hence 
the reason for the factor of 1/2 to avoid doublc-coimting of the contributions from the 
diagonal terms — is clear. We note that for s e N, s > 2 the numerical value of the sums 
is independent of any manner of grouping terms. Even so, the form of the integrands in our 
integral representations reflect this choice. We have found that the convention (14) yields the 
most symmetric expressions in appearance (a different grouping for integer s was employed in 
[7]). 

Next, we recognize that in sums of the form (14) we may substitute the appropriate plane- 
wave expansion (12) to represent the poles contained in /. This substitution will transform 
the quadrant sums into geometric series. Concerning the later, we derive the following lemma. 

Lemma 2. For any p,q G C and K gN, the following is an identity 



K i , . -1 , „ , N , /, , X r/ -l^K+1 /„n^K+l 



V S^e vW - 1 p(g-^+2 + g) _ 1 ( 1 + q 
fer.e, 2(l-?>(g + «-i)+p^) 2[l-q 



1 — pq ^ 1 — pq 



(16) 



Proof: 

The formula follows from iteration of the usual single variable geometric sum, and algebra. 

We record the following corollary for reference as the expressions appear many times in the 
formulas below. 

Corollary 3. We have the following specializations of lemma 2 

K m 

J2 E e™„(e-^)'"(e-^-)" = 2e-Vi(T, A) - 2e-^(^+i)/}^)(T, A) 

m=l n=— m 

and 

K n 

E E en-K^)"(e'^)'" = 2e-V2(r, A) - 2e-^(^+i)/i^)(r, A) 



n=l 171=— n 
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where the functions /i, f[^\f2, 12^' are: 



/i(r,A) 



h{T,x) 



cc)slr(rA/2) 



1 - 2e-^ cosh(rA) + e" 



1 /I + e- 



-2A 



-\t{K+1) 



1 



4 Vl-e-^ 

cos^(A/2) 
1 - 2e'-^^ cos(A) + e2i-r^ 



-a(i-t) 1 _ e-A(i+T) 



1 / 1 + c'^ 
4 I 1 - e'-^ 



-iA(_ff+l) 



1 _ QiX(,T-l) I _ eiA(r+l) 



(17) 

(18) 
(19) 

(20) 



We assume r is in the fundamental region (10) and make several observations. First, we 
note that all of the functions given by (17)-(20) have double poles at the origin, A = 0. Since 
Sy(T) > 0, neither /i(t. A) nor /2(t, A) have other poles for A > 0. For r strictly imaginary, 
the denominator 1 — exp(— rA) of /{^'(t, A) will have isolated simple zeros. However, these 
are balanced by simple zeros of the difference of bracketed terms in (18). Thus the product, in 
other words ,f[^\ has no other singularities for A > 0. A similar argument shows that f2^^ is 
also finite for A > 0. Finally, one may show the bounds 

|e-^/i(T,A)| < Cie-^(i-l^(^)l), and 
|e'^V2(T,A)| < C2e'^^ 

for large A. As |5i(r)| < 1/2 and 3(t) > 0, both quantities arc exponentially decreasing in 
A. Similar reasoning shows that e~^^^'^^'> fi^\T,\) and e''^'^*^''^"'"^)/2^'*(r. A) are exponentially 
decreasing in A and K. 



3 Eisenstein series 

As mentioned previously, the summation for 3?(s) > 2 is absolutely convergent. We begin 
by proving our first integral representation for this case in theorem 4. As a corollary, by 
restricting s = n, n > 3 wc obtain integral representations of _E„. Further inspection of the 
integral representation demonstrates the existence of £"2- Elaborating on theorem 4, we derive 
an alternative representation for Es as a contour integral. As a consequence of this second 
representation, wc prove that Es admits an analytic continuation in s as an entire function. 
Returning to the analysis of E2 , we consider a more general limiting procedure and define £^2"'' 
as the limit over increasing "lattice rectangles" with a fixed aspect ratio defined by a. We 
derive a closed form expression which, when added to E2 = E'ip, gives eI^\ As a corollary, 
we derive the relationship between E2 and the sum E2 as defined by Eisenstein. 

3.1 Integral representations 

For the sums Es defined by (5) we prove 
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Theorem 4. Given a lattice A(t) with raiio t chosen from the fundamental region (10), we 
have the following integral representation for Eg, 3?(s) > 2 

E,{t) = cos ^ (e-'^^/^e-Vi(r, A) + e^VsCr, A)) dA (21) 

where /i(t, A) and /2(t, A) are given by (17) and (19). 



Proof: 

We recognize that, due to the placement of branch cut (6) and the symmetric summation 
conventions (14), we have the following relations between sums over and _|_^ 



,,,s ^ ,,,s 



Therefore we may consider the positive quadrants only and scale the results by an exponential 
factor. Turning to the quadrant A^,^, in place of the isolated singularity of degree s, we 
substitute the appropriate plane- wave expression from (12) to obtain 



K ni 



y 7 — ^ = {i+e~'ny y 



(m + nrY ^ (m + nrV 

'^^^(± .) n=-m ^ ' 

K 



(e-Vi(T,A) -e-^(^+i)/f )(T,A)) dA 



2(1 + e- 



T{s) 

where the last line follows from corollary 3. From the statements following this same corol- 
lary, we observe that the two integrands are singular at A = 0, and are otherwise finite and 
exponentially decreasing in K and A € R"'". In addition, as 3?(s) > 2, the singularity at the 
origin is absolutely integrable. Therefore, one may take the large K limit inside the integral 
and compute 

lim y = cos (Is) r A"-^e-'"'^/2e-^/i(r, A)dA. 



By a similar analysis we prove that 



lim y- = cos f ^s) / A"-ie'^^/2(r,A)dA. 

K^oo ^ (m + nr)^ \2 J T{s) Jo "'^^ ' ' 
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Adding these two contributions gives the theorem, f 



A few comments are in order. First, for s = 2j + l the term cos (§s) = 0. This is equivalent 
to the well-known fact that, by symmetry, the odd Eisenstein series are zero. Furthermore, by 
inspection of (21), we see that the absolute convergence of the Eisenstein series i?s,5ft(s) >2 
manifests itself in the behavior of the integrand of (21) near the origin; the factor A"^^ balances 
the double poles of /i and /2 so as to ensure the product is integrable at A = 0. More careful 
analysis reveals that the formula (21) is finite even for the conditionally convergent case £2- 
The Laurent expansions of the integrands about the origin are 

cosh^(rA/2) ^ 1 _ (1-2t^)A , ^ 3n 

A "0 1 - 2e-^ cosh(rA) + c-2A A(l - r^) 12(1 _ r2) ^ 

cos2(rA/2) ^ 1 _ (2-r^)A . 3^ 

A^o 1 - 2ei^A cos(tA) + e'2^A A(l - r2) 12(1 - t2) ^ 

At s = 2, where the exponential factor e~'*'^/2 — ^j^g expansions above are to be subtracted. 
Hence the integrand of (21) is finite at the origin even in this case. By a similar analysis one 
may show that the if-dependent terms also cancel at the origin. We have proved: 

Corollary 5. The summMion (5) converges in the conditionally convergent case s = 2, and 

its value, E2, is given by the integral (21). 

In fact, a great deal more may be said. The function 

^/ ^ io^/o , cosh^(T2;/2) • ^ cos2(z/2) 

Fis z) = e ' e~ - - he — - f22) 

^ ' ' 1 - 2e-^ cosh(r2;) + e-2^ 1 - 2e''^^ cos(^) + e2i-^^ ^ ' 

appearing as a factor in the integrand (21) has a singularity at the origin z = 0, and additional 
simple poles in the complex plane at the points 



27ri 27r 
z & P = { ±- m, ±- 



l±r ' 1±T 



m, n G N 



We denote the minimum magnitude of a\\ z & P by p. Next, define the contour C which begins 
at 00 + iy, y > 0; runs parallel to real axis until it intersects the circle centered at the origin 
with radius r, where y < r < p; follows this circle counterclockwise around the origin; and runs 
back out to 00 — iy, parallel to the real axis. We assume that y > is small enough such that 
C encloses only the pole at z = 0. Finally, for the function z^~^,s ^ Z, situate the branch cut 
along the positive real axis such that (A e R) 

lim(A + i2/)*-^ = X'-'^ 

lim (A - iyy-^ = e2'^i(«-i) A«-i ^^^^ 
With these preliminaries established we prove the following theorem. 

Theorem 6. Given a lattice A(t) with ratio r chosen from the fundamental region (10), we 
have the following contour integral representation for Eg, 5i(s) > 2 

r(i -^c 

i-K 



Eg = 2 cos I 
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where F{s,z) is given by (22). 



Proof: 

Consider the contour integral 



L 



z'-^F{s,z)dz. 



As the integrand is analytic except for the singularity at the origin, we may deform the contour 
without altering the value of the integral. Specifically, shrink the radius of the circle, r — > 0, 
and take the limit as y ^ for the two components running parallel to the real axis. For 
3fi(s) = 2 + e, we estimate the contribution from the circular arc 



lim 



/ z'-^F{s,z)dz <M\im ( r'" V-^l = 0. 

J\z\=r '•^0J|2|=r 



Turning to the components parallel to the real axis, applying the specification of the branch 
cut (23) we compute 

reT+\y poo 

lim/ z'-'^F{s,z)dz = - X'-'^F{s,\)dX 

y^'^ Joo+iy JQ 

roo—'iy pco 

lim/ z'-'^F{s,z)dz = e'^'''^'-'^^ X'-'^F{s,X)dX. 

V^OJe^-iy Jo 

We recognize the integrals of theorem 4, make use of the identity 

TT 



r(s)r(i -s) = -: 



sin(7rs) 



and obtain 



/ z'-'^F{s,z)dz = (e2'^'('*-i) - 1) f X'-^F{s,X)dX 
Jc Jo 



2ie'^''' sin(7rs)r(s) ~ 

7 \ -^8 

4 COS (|s) 
in 



2cos(fs) e-i«'^r(l-s) 
The result (24) follows from algebra, f 

Several corollaries follow from theorem 6. We note here only the most immediate 

Corollary 7. The sums Es, defined by (5) for 5R(s) > 2, admit an analytic continuaMon to 
s € C as an entire function. This continuation is given by the contour integral representation 
(24). 
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Proof: 

The contour integral appearing in (24) defines an analytic function of s which is never singular. 
The same may be said for the cosine factor. Thus the only possible singularities would arise 
from the factor r(l — s) which has simple poles for s G N. From the definition (5) we know 
that Eg is finite for 3?(s) > 2 (the apparent singularities in (24) in this case are balanced by 
zeros of the cosine term, the contour integral, or both.) Therefore, we need only verify the 
finite existence of Ei and i?2- We have argued above that E2 is finite. Finally, for s = n = 1, 
the pole in the Gamma function is balanced by the simple zero of the cosine factor, f 

By inspection of (24), we find that Eg has simple zeros for s = 1 — 2j, j > 1. Thus 
the continuation respects the symmetric limiting process (5). We anticipate further results 
concerning the evaluation of the contour integral (24) via residue methods. However, we have 
not carried out this analysis at the time of writing this paper. Finally, as mentioned previously, 
the exact form of the integrands (21) and (24) reflect our summation convention with respect to 
grouping of summands and placement of the branch-cut. Regarding the latter, similar formulas 
arise if the branch-cut is situated along any of the lattice diagonals; the effect is to redistribute 
factors of exp(i7rs/2) between the two functions /i and /2. There are, perhaps, additional 
treatments of the branch-cut that could yield relatively simple expressions. However, a simple 
integral expression valid for placement of the cut along an arbitrary ray in the complex plane 
appears to be intractable. 



3.2 Aspect ratio correction 

Although we have demonstrated that E2 is finite, until this point its relation to Eisenstein's 
definition of this series is unclear. We make this explicit in the present section. 

We wish to formalize summation over lattice rectangles with a fixed aspect ratio. Given 
a e (0, 00) we define the sum 

= hm y y (25) 

\m\<yaK\ \n\<K ^ ' 

where [a;J denotes greatest integer less than or equal to x. Again, although expected, the 
existence of the limit (25) is not a priori guaranteed but will follow in the course of our 
analysis. 

Clearly, e'^'' = E2 defined in (5). More generally, we write 

=^2 + A(a,T) (26) 

where the value of E2 may be computed via the integral expression (21). Analysis similar to 
that used to prove (21) may be employed to compute a closed form expression for A(a, r). 

Theorem 8. For a fixed aspect ratio determined by a & (0,oo), the limit E^'^ specified in (25) 
exists. Furthermore, when written in the form (26), the aspect ratio dependence is given by 

A{a,T) = — — ^arctan(iT) — arctan (^'^) 
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Proof: 

Assume a > 1. We write the limit (25) 



laK] n=K 

A(a,r) = 2 1im E 



m=K+ln=-K ^ ^ 

The contribution from the sum over lattice points — [a ■ K \ < m < —K — 1 is accounted for by 
the factor of two multiplying the sum in the final line. We represent the poles using the 
plane- wave expansion (12). 

yaK\ n=K 



A(a,T) = 2 lim ^ ^ / Ae-^(™+"^)dA 

m=K+l n=-K •' ^ 



= 2 lim / 



qXtK _ ^-\t(K+1) \ /q-\(K+1) _ ^-\\_aK\ 



Kt _ -Xt(1+1/K) \ / -A(l+l/_ff) _ ^-\yaK\/K ' 

2 lim / A ( ^ . ^ , ( : %- ) -^^dA 



2 r(e^^-e--)(e--e--)^^_ ^^^^ 



T Jo A 

Where the argument which justifies taking the large K limit inside the integral runs along the 
same lines as in the proof of theorem 4. This last integral (28) may be evaluated in closed form 
using the formula 

J e"'^^ sm{Sx)-dx = arctan l^-j , (29) 

which holds for 3?(/?) > |S((5)| (sec [4], 3.944.5). Taking care to write (28) as the difference of 
two integrals of the form (29), and performing algebra gives the expression for A(a, r), a > 1 
in (27). 

For a < 1 the lattice rectangle is such that the longer side is in the r-direction. As written, 
equations (25) and (26) suggest that this rectangle is inscribed in a lattice square of size K, 
and to compute A(a, r), one should subtract the extra contributions exterior to the rectangle 
but interior to the square. The problem with this approach is that, for arbitrary r and a, one 
would have to keep track of the quadrants in which these points lie. In lieu of this, for a < 1 
we rescale the limits in (25) 



E. 



f)= lim y y 



{yify -(- 77,7"^^ 
\m\<K\n\<\KI<x\ ^ ' 

Informally, this has the "effect" of inscribing the square in the rectangle and motivates com- 
puting the contributions from the points in the difference using the plane-wave formulas ap- 
propriate for A^ _|_^ . Arguing as above and using the integral identity (29) we compute 

A(a,r)=-2 1im EE/ Ae'^^""^"^-'dA 

n=K+l m=-K ■'^ 
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qA 



r 



A 



= — — ^arctan ^"'^^ ~ arctan • (30) 

Although perhaps not obvious at first glance, the formula (30) is the same as the formula 
for A(a, r) in (27). Using standard trig identities we have 

lim ^arctan(z) — arctan ^ = lim arctan^ — — ^ = ±-^. (31) 

Furthermore, for z = ir or z = it /a with t satisfying (10) and a > 0, we find that we should 
choose the minus sign in (31). Using this identity and algebra we observe that (30) is equal to 
(27), thus (27) holds for all a > 0. f 

We use this theorem to find the connection between our summation and Eisenstein's. Start- 
ing from Eisenstein's summation convention we obtain 

N M ^ 

Eo = lim lim -- 

jv^oo M^oo (m + nr)^ 

n=-N m=—M ^ ' 

N m=—[a-N\ ^ 

= lim ^ lim -- 

N^oo ' a—^oo ^ — ' [171 + riT) 
n=-N m=-la-N\ ^ ' 

N m=[a-JVJ _^ 

= lim lim -^r 

a^oo Af^oo ' ' [m + tit) 

n=-N m^-la-N} ^ ' 

4i 

= £2 arctan(iT). (32) 

r 

Standard estimates justify commuting the N and a limits between the second and third lines, 
and wc used (27) to compute this limit. As r 7^ 0, (32) shows that, for finite t, the value of our 
sum is always different from Eisenstein's. In the limit r ioo, |3?(r)| < 1/2, both summations 
are equal and presumably converge to 2(^^(2) = 7r^/3. 

In a similar vein, wc compute the difference between taking Eisenstein's limit and "its 
reverse" . Arguing as above we have that 

M ^ \ ~ 

lim lim y^ = £2 + hm A(q:,t) 

m=-M n=-N ^ ' 

= E2-^ (^arctan (^^^ + (33) 
Taking the difference between (32) and (33) and we obtain 

/ N M M N 

lim > lim > — lim > lim > 



1 _ 27ri 

N—i-ca ' M—i-ca ' M—i-ca ^-^ N—i-ca ' / (m + TIt)^ T 
I n=—N m= — M m= — M n=—N/ ^ ' 
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For a different proof of this fact see Walker, [11]. 

Finally, in the case of the square lattice, r = i, we observe that /i(i. A) = /2(i, A). Collecting 
factors in (21) and performing algebra we find that 

Q 1*00 

^,(i) = e-'-/4 cos cos (|.) A^-^e- A)dA. 

Corresponding to the added 7r/2-symmetry of the square lattice, the product of cosines causes 

the sum to vanish for s = n > 3, n not a multiple of four. Similarly, as the representation (21) 
is valid for n = 2, we conclude £'2(1) = 0. Substituting this value into (26), and taking the 
large a (small a) limits we obtain 

JV M ^ 

lim ^ lim ^ — r- = tt 

n=—N m= — M ^ ' 

M N ^ 

lim >^ lim >^ 



M^oo Jv-»oo ^-^ (m + in)^ ' 

m=-M n=-N ^ ' 

well-known identities in the fast multipole community (see, for example, [5]). 



4 Weierstrass elliptic functions 

Our derivations of integral formulas for Weierstrass' elliptic functions proceed in much the same 
manner as above. As with the Eisenstein series, wc begin with the definitions of the functions 
as sums over the lattice (7) and (8). Next, we group terms in the sum as in (14), substitute 
the appropriate plane- wave expansions from lemma 1, and sum the resulting geometric series 
using the identities of corollary 3. 

As a preliminary note, the integral representations for p(2, r) and <^(z,t) which we derive 
in theorem 9 are not valid for all 2: € C, but rather have a finite domain of validity. This is 
a consequence of the way in which we group terms. More precisely, we require that z ^ D{t) 
defined by 

D{t) = {z\ 3fi(-l ± z ± r) < 0, 9(t ± 2) > 0}. (34) 

As T is in the region (10), one may verify that D{t) is an open set containing the origin. 
However, wc point out that D{t) may not contain the fundamental period parallelogram of the 
lattice. For example, the standard hexagonal lattice has generators (1, r) = (1, 1/2 + 1-^3/2). 
Thus a corner of Ag is given by the point zq = 1/2 + r/2 = 3/4 -|- i-\/3/4. However, JR(— 1 -|- 
-20 + t) = 1/4 > 0, violating the first inequality in (34). 
With this aside, we prove the following. 

Theorem 9. Assume A = A(/x, v) is an arbitrary complex lattice with generators chosen such 
that T = ly/fi is in the fundamental region (10), and that the complex number z is in the 
domain D{t) defined by the inequalities (34)- We have the following integral expressions for 
the inhomogeneous elliptic functions p{z,t), and ({z,t): 

1 f°° 
p{z,t) = — +8 a 
z Jo 



' sinh^ 



zX 



/i(A,r)+e'^^sin- 



zX 



/2(A,t) 



dA, 



(35) 
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and, 



1 f°° 

({z,t) = -+4 / [e-^(^A-sinh(2A))/i(A,T) -e'^^(^A-sin(2A))/2(A,T)] dA, (36) 
^ Jo 

where the functions /i,/2 are defined by (17) and (19). We may evaluate the homogeneous 
functions, p{x\ijl,u) and Q{x\ii,v), via the appropriate scaling relations (9) and (35)-(36) with 
the proviso that x/ ii = z & D{t). 



Proof: 

In computing the integral representation for the p function we will group terms of the sum (7) 
as in (14). As in the computation of the Eisenstein sums, we wish to combine the contributions 
from the sums over . We compute 

5 \{z-ujy ~Z^J^^ ^ ^""^ \iz + m + nTY ~ (m + nrf 

^ \(z-ujY ~ ^) " ^ ((-^ + m + nr)2 ~ (m + utY ) ' 

<^gA^ m=ln=-m ' ^ ' ' 

Therefore the contributions from both quadrants may be expressed as a single sum over A^ 
of a modified summand. Furthermore, under the assumption z & D, all of the poles in this 
sum may be expressed using the 3?(a;) > plane- wave expansion from (12). 

^ ( 1 1 A V V / 1 2__ 1 

^ \iz-ivy to^J ^^'""V(-^ + m + nT)2 (m + nr)^^ (z + m + riT)^ 

/.oo /Km \ 

= I A(e^^-2 + e-^^) ^ ^ e„„e-^('"+"-) dA 

yrn— 1 71— —m / 

= 8 Asinh^ (e-Vi(r, A) - e-^(^+i)/f ^(r, A)) dA. (37) 

Arguing as before, wc find the large K limit of the i^-dependent term to be zero. We compute 
the contribution from the terms in the quadrants A^ _|_^ in an analogous manner. Adding this 
result to (37) gives (35). 

The derivation of the expression for the C function is similar. In brief, the sum over the 
quadrants A^ may again be expressed as a sum over the single quadrant A^ in which we 
substitute the appropriate plane-wave expansion. Thus, 

V-/1 ^ z\ ^ ^ / 1 2z 

1^ 7 T ^ ^^1 = 1^ V Cmn ■ ■ h ■ 

\K \ \ / / TO=1 7l= — m ^ ^ ' 



(±,») 



^00 /Km \ 

= / (e^^ + 2^A - e-^^) ^ ^ e„„e-^("+"^) dA 

\m=ln=-m / 
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POO . . 

= {z\ - sinh(0A)) (e-Vi(A,r) - e-^^^+^VrV^r)) m) 

As before the -ftT-dependent term goes to zero in the hmit. The analogous sums over _|_^ 
give the other half of the expression (36). f 

Remarks: 

1. As an alternative to the above derivation of theorem 9, we recall that the Eisenstein 
series appear as coefficients in the Laurent expansion for Weierstrass' p function 



wGA\{0} 

= \+ y y(n -1)^-2 J- 

wGA\{0} "=3 

oo 

n=3 

Substituting the integral representations (21) for the E^^ the Taylor series may be 
summed explicitly inside the integrand. The formula (35) above follows after algebraic 
simplification. Furthermore, the expression (36) for the C function follows from anti- 
differentiation of (35). 

2. As with the series Ei, and its dependence on aspect ratios derived in theorem 8, the 
slowly decaying terms of the sums defining p and C manifest themselves at the origin 
in the integral representations (equations (7), (8) and (35), (36) respectively). In the 
integral representations, we observe that Weierstrass' "correction" terms are arranged 
in such a way as to create third order zeros at A = 0, which appropriately balance the 
second order poles from /i and /2. 



5 Conclusion 

We conclude with a brief discussion of our results in relation to previous research in this field. 
To the best of our knowledge, there is no analog to the integral expressions for the p and 
C functions (35) and (36). The possibility of developing numerical routines for evaluation 
of these functions based on these representations deserves further study. We observe that 
the integrands arc not extremely oscillatory, and decay exponentially. Thus A'^-point Gauss- 
Laguerre quadrature rules will converge rapidly in A. As one drawback, there is the perhaps 
awkward domain of validity in z. However, it may be that symmetries of the p and C functions 
imply that it is sufficient to evaluate them over domains that are smaller than the fundamental 
period parallelogram. Furthermore, at least for the p function, there exists the following 
closed-form Fourier-like expansion [10] 

p(z,r) = -2 - + y + 87r2 V , . cos(27rnz) . (39) 
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Both summands in (39) arc exponentially decreasing and the sums converge rapidly stiff 
competition from a numerical perspective. Nevertheless, we have not fully explored the relative 
merits of this approach over the plane- wave representation (35). In addition, the integral 
representations may have further analytic implications. 

Turning to the representations for the Eisenstein series, the existence of E2 (corollary 5) 
is not unexpected. In addition to the original finiteness proofs given by Eisenstein, many 
years prior to this present work. Walker derived the remarkable formula for the conditionally 
convergent series (see [10]) 

lim > 7 ^ = : 47ri , , , 

where the Dedekind 77(T)-function with 9(r) > is defined by 

00 

r^(r) = e"'"/i2 ]J (1 _ e^^^'^"). 

We also note that a different treatment, initiated by Heckc, has become a standard approach 
to resolving convergence and transformation properties of E2 [9] . 

As indicated by theorem 6, our expressions are quite general, and have broad implications. 
We note that Riemann demonstrated both the fimctional equation satisfied by CRi-^)i s-^id the 
evaluation of C,R{—2n + 1) (and, via the functional equation, ^ij(2n)) in terms of Bernoulli 
numbers using the "version" of theorem 6 appropriate for his^zeta function. Similarly, we 
anticipate that a residue argument will give the evaluation of En ~ En in terms of multiple 
Bernoulli numbers. For an alternative treatment of Eisenstein series for negative even integers 
using Hecke convergence factors see the recent work of Pribitkin [8]. The functional equation 
satisfied by the continuation of Es is more elusive. We are currently pursuing this and hope to 
report our results in the future. 

Finally, there is a possibility that representations of the form (21) may exist for certain 
Dirichlet series 



G{s,x)= E 

i^eA\{o} 



x(w) 



x(a;) = exp(i(TO/(ia + ny(3)), 

for a, /? G R. (These arc called "Kroneckcr series" in [12], Chapter VIII.) A detailed discussion 
of the convergence of these series is given in [2] . We note that Laplace-Mellin techniques have 
been employed frequently in the analysis of such series. The approach up until now has been 
to think of 

'' ' — [m/i + nv]"^ — Q{'m, n) 



-^m.n I 



as defining a positive definite quadratic form taking m and n as arguments. Treating this form 
as "indivisible", one may use the 3?(a;) > plane-wave formula in lemma 1 and obtain the 
integral representation 

V 2^ = J_/'"a«-i V Y(a;)e-^'3(™'")dA. 
^ |a;|2« r(s) Jo ^ ' 
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The analysis then proceeds via 9 functions. 

Our app roach w ould be different. "Plane-wave-like" representations exist for the function 
f{x, y) = 1/^^"+^. Formally, one may take the true plane-wave expressions for f{x, y, z) = 
\J + ip" + derived in [3], and set z = 0. The result is a 2D integral — as opposed to the 
Eisenstein case analyzed above where one complex dimension (two real) collapses into a ID 
integral. However, the critical element of this representation is that the exponential function 
in the "plane-wave" representation is linear in m and n. Again, as a consequence, we observe 
that the summation under the integrand becomes a 2D geometric series. We are considering 
this as a possible direction for future research. 

Acknowledgments. The authors thank Peter Walker at the University of Sharjah, United 
Arab Emirates, and Brad Alpert at NIST/Boulder for frequent helpful discussions during the 
writing of this paper. 
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Figure 1: Partition of into subregions. The generators (1, r) are shown in red. The central 
dotted region is the boundary of the fundamental domain. The dashed lines show the divisions 
into K¥, , The solid black line is the branch cut. 
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